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Due to its large mass, the top quark has the largest de- 
cay width and hence the shortest lifetime of the quarks 
in the standard model (SM). The total width of the 
top quark at the leading order is dependent on the top 
quark mass m t and the Fermi coupling constant Gf'- 
r° = Girm|/(87rv2)- Higher order effects include 
the introduction of finite W boson and b quark masses 
(Mw,Trib), non-zero off-diagonal elements of the quark- 
mixing matrix, and higher order corrections in the strong 
coupling constant a s . Neglecting terms of order m\jm\ , 
Qj, and (a s /7r) M^/m 2 , the width predicted in the SM 
at next-to-leading-order is: 



TABLE I: Event selection requirements for the 1-tag and 2- 
tag event samples. 



r t = r?(i 



^t) 2 (i 



2^1 



)[1 



2a s 2tt 2 5 
IWr IT ~ 2"' 



The total width of the top quark is calculated to a preci- 
sion of about 1% in the SM; it is approximately 1.5 GeV 
for m t = 175 GeV/c 2 [HQ. 

A deviation from the SM prediction could indicate a 
significant contribution from top quark decays to non- 
SM particles such as t — > bH + (where H + is the charged 
Higgs boson in the supersymmetric model), or from rare 
SM processes such as t ^ dW + and t — > sW + . Al- 
though such scenarios have not been observed experimen- 
tally [H, 0, H, S 0] , a general way to rule out the presence 
of a large top quark decay rate to non-SM channels, in- 
cluding those with non-detectable final states, is through 
experimental constraints on T t . To date, there have been 
no direct experimental measurements of the total width 
of the top quark. 

The data set for the analysis presented in this paper 
is collected by the CDF II detector, a multipurpose par- 
ticle detector for pp collisions at y/s — 1.96 TeV at the 
Fermilab Tevatron. A charged particle tracking system 
immersed in a magnetic field consists of a silicon mi- 
crostrip tracker and a drift chamber. Electromagnetic 
and hadronic calorimeters surround the tracking system 
and measure particle energies. Drift chambers and scin- 
tillators located outside the calorimeters detect muons. 
The detector is described in detail elsewhere Q. 

We employ a cylindrical coordinate system where 9 
and <f> are the polar and azimuthal angle, respectively, 
with respect to the proton beam. Transverse energy and 
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Event selection category 


1-tag 




2-tag 


E T (GeV) or p» (GeV/c) 




> 20 




St (GeV) 




> 20 




Jets 1-3 E T (GeV) 




> 15 




Jet 4 E T (GeV) 


> 15 




> 8 


Number of b tags 


1 




> 2 











momentum are Et — EsinQ and px — \p\sinO, respec- 
tively, where E and p are energy and momentum. Miss- 
ing transverse energy, St, is defined as the magnitude of 
the vector — 'SiE T n- 1 where E T is the magnitude of trans- 
verse energy contained in each calorimeter tower i, and 
rii is the unit vector from the collision point to the tower 
in the transverse plane. 

Top quarks are produced primarily by strong inter- 
action in top-antitop (tt) pairs at the Tevatron. Top 
quarks decay almost exclusively to a W boson and a b 
quark through the weak interaction in the SM. We iden- 
tify candidate tt events in the "lepton + jets" channel, 
where one W boson decays to an electron or a muon, 
and a neutrino, while the other W boson decays to a 
quark-antiquark pair. We select events consistent with 
this topology, requiring a high-p-r electron or muon can- 
didate, missing transverse energy denoting the presence 
of a neutrino, and four jets. Jets are reconstructed using 
a cone algorithm with radius AR = \J Arj 2 + A<j) 2 =0.4. 
At least one jet with Et > 15 GeV must be identified 
as a b quark candidate through the presence of a dis- 
placed vertex within the jet cone arising from the decay 
of a long-lived bottom hadron (fe-tag). The event selec- 
tion criteria are listed in Table HI Detailed information 
on event selection is available elsewhere 0, [l(| ■ 

We divide the candidate events into two exclusive 
classes: one (1-tag) containing events with one 6-tagged 
jet among the leading four, and another (2-tag) with two 
or more such jets. Separating these subsamples results in 
a more efficient use of statistical information due to their 
different reconstructed top mass resolution and signal- 
to-background ratios. Background events are expected 
primarily from W production in association with jets 
(W + jets), multijet processes where a jet is misidentificd 
as a charged lepton and St results from energy mismea- 
surement of the jets (non-14 7 ), and small contributions 
from electroweak backgrounds (EWK) composed of sin- 
gle top quark and diboson (WW, WZ, ZZ) production. 
Table [TT] summarizes the expected sample compositions 
that are obtained by scaling to 955 pb _1 from a previous 
tt analysis with 318 pb _1 [11[. A detailed description of 
the background estimation is given in Ref. [9j. 

After event selection, the analysis proceeds in three 
steps. First, we reconstruct a top quark mass m r t cco from 
each event using a kinematic fitter. The width of the 
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TABLE II: The sources and expected numbers of background 
events, and the number of events observed for the 1-tag and 
2-tag event samples in our 955 pb _1 dataset. 



E 

o 

f= 0.02 
5 2-tag signal shapes 



Background source 


1-tag 


2 


-ta 


g 


W + jets 


21.4 ± 4.9 


4.6 


± 


1.5 


non-W 


7.0 ± 1.5 


0.9 


± 


0.2 


EWK 


2.8 ± 0.6 


0.17 


± 


0.04 


Total background events 


31.2 ± 7.0 


5.7 


± 


1.8 


Observed events 


171 




82 





reconstructed mass distribution for the selected events is 
sensitive to T t . The second step is a likelihood fit of the 
reconstructed mass distributions using simulated signal 
and background distributions that yields an estimator 
of T t . Finally, we use a frequentist prescription (with 
Bayesian treatment of systematic uncertainties 121]) to 
determine a 95% C.L. upper limit on T t in the physically 
allowed region. 

We perform a x 2 minimization to fit the momenta of 



the tt decay products and determine m\ cco for each event 
using the four leading jets. We assume that the final 
state arises from the decay of a tt pair into W bosons 
and b quarks. To resolve the ambiguity arising from the 
different ways of assigning the jets to the four quarks, we 
require that 6-tagged jets are assigned to b quarks and 
select the assignment with the lowest \ 2 ■ This kinematic 
fitter is used in other CDF analyses and is described in 
detail in Ref. [U|. In the x 2 fit, both sets of W decay 
daughters are constrained to have the invariant mass of 
the W boson, and both Wb states are constrained to have 
the same mass. Although the top and antitop quark will 
likely be produced with different masses, we confirmed 
that there is no significant difference in the sensitivity 
resulting from the m t = m t - condition even for a large 
value of T t . 

To distinguish between different values of T t , we com- 
pare the m T t ec ° distribution from our data to a series of 
samples created using the PYTHIA 6.216 event genera- 
tor [13| and a full detector simulation. We use samples 
with m t = 175 GeV/c 2 and various values of T t between 
0.001 GeV and 100 GeV. Although pythia does not fully 
account for quantum interference with irreducible back- 
ground diagrams and off-shell effects, for T t < 30 GeV 
these effects are small, and the existing description is 
expected to be adequate [14j. W + jets background 
events are generated using ALPGEN 1.3 [15f . with par- 
ton showering and fragmentation in herwig 6.504 16| . 
An unbinned extended maximum likelihood fit 17] is per- 
formed using parameterized signal and background m\ ec ° 
templates. As an example, Fig. Q] shows templates for 
the signal m\ eco distribution in the 2-tag subsample at 
three values of T t . We parameterize the m\ eco distribu- 
tions as a function of T t . Small shifts in the mean of 
the templates are induced by the interplay between the 
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FIG. 1: The parameterized signal m\ cc ° distributions and 
goodness of the parameterization are shown for the 2-tag sub- 
sample at three different values of Tt. The parameterization 
is determined in a global fit to all the 2-tag templates. The 
parameterized signal mJ eco distributions for 1-tag are similar. 



top mass Breit-Wigner distribution and the parton distri- 
bution functions that preferentially produce events with 
low quark masses. In the likelihood fit, we constrain the 
background templates in the 2-tag and 1-tag samples to 
the levels given in Table HI The best fit value if* is the 
width which maximizes the likelihood. 

We allow negative Tf * values that represent m\ ec ° dis- 
tributions narrower than the nominal due to statistical 
fluctuations. The expected m\ eco distribution for a neg- 
ative is derived from an extrapolation of the param- 
eterization to the unphysical region. The reconstructed 
top mass distributions from the data and the results of 
the likelihood fit are shown in Fig. [5] The data are con- 
sistent with the fitted curves with the preferred value of 
rf* = —4.8 GeV. For the data sample used in this anal- 
ysis, we expect to measure a negative rf* about 40% of 
the time and rf* less than -4.8 GeV about 25% of the 
time if the total width is 1.5 GeV. The limit on the true 
value of r t , however, will be constrained to the physical 
region. 

To set a limit on the true value of T t , we employ 
the Neyman construction [l8| to ensure a coverage of 
at least 95%. The likelihood-ratio ordering principle due 
to Feldman and Cousins (l9j provides a smooth transi- 
tion from one-sided to two-sided limits and usually guar- 
antees a non-empty interval. We derive the confidence 
belts from ensembles of simulated experiments in which 
signal events are selected from the simulated samples gen- 
erated with different values of IV The rf* distribution 
from each such ensemble is convoluted with a shape that 
represents the effects of systematic uncertainties as de- 
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FIG. 2: The m\ ac ° distribution is shown for 1-tag and 2-tag 
data samples overlaid with the signal and background distri- 
butions from the combined fit. A rf* likelihood scan is shown 
in the inset; for the shaded non-physical region (Tf* < 0), the 
likelihood uses m\ ec ° distributions extrapolated from the pa- 
rameterization to shapes with Ff 4 > 0. 



scribed below. 

Since the top quark mass reconstruction is dominated 
by the measurement of jet energy, and since our fit is 
largely determined by the peak and the width of the m r t eco 
distribution, the uncertainties on the jet energy scale and 
the jet energy resolution are the dominant uncertainties 
in the top width measurement. The uncertainties on the 
jet energy scale calibration are extensively studied using 
a combination of simulated and data control samples [20| 
and amount to about 3% of the measured jet energy for 
jets in the ti sample. The effect on the rf * distribution 
is non- linear, with larger T^' being more likely for larger 
jet energy scale shifts in either direction. This is because 
the only degree of freedom in the templates is the width, 
and a signal template with a larger T t accommodates the 
events with the shifted peak. 

We select events with one jet and one high-pr pho- 
ton and compare their energies to study the jet energy 
resolution. Data and PYTHIA events show similar jet 
resolution of 15% - 10% for jet transverse energies be- 
tween 20 GeV and 200 GeV, respectively. Taking into 
account statistical uncertainty of the data, we define a 
PT-dependent systematic uncertainty on jet resolution of 
10% - 4% to cover the difference. Then, we add Gaussian 
smearing with corresponding uncertainty to each jet in 



signal Monte Carlo events. We also study smaller system- 
atic uncertainties in T^' related to the background m\ ec ° 
shape, Monte Carlo statistics, the Monte Carlo genera- 
tor, the parton distribution functions, and other signal 



modeling effects [111 ]. The combined convolution shape, 
accounting for all systematic uncertainties, has a shift of 
-0.4 GeV and an RMS of 4.6 GeV. This can be com- 
pared to the RMS from statistical effects only: between 
6.6 GeV and 10.1 GeV for simulated experiment ensem- 
bles using T t of 1.5 - 30 GeV. Figure 02 shows the 95% 
confidence belts after including systematic uncertainties. 
The fitted value from data, Tf * = —4.8 GeV, corresponds 
to a limit of T t < 13.1 GeV at 95% C.L. 

Our measurement assumes a fixed value for the top 
quark mass of 175 GeV/c 2 . The one-dimensional like- 
lihood is sensitive to this assumption in the same way 
as described above for jet energy scale uncertainties. In 
particular, if the top quark mass is consistent with the 
current world average of 171.2±2.1 GeV/c 2 2l|, the con- 



fidence belts would shift to higher 1^ values, resulting 
in an upper limit of T t lower than what we quoted in this 
paper. 

In summary, using 253 top-antitop pair candidate 
events we present the first direct experimental up- 
per limit on the total decay width of the top quark, 
T t < 13.1 GeV at 95% C.L. for m t = 175 GeV/c 2 . 
This corresponds to a limit on the top quark lifetime of 
r t > 5 x 10 -26 s. This measurement is statistically lim- 
ited and its dominant systematic uncertainties are likely 
reducible with statistics. The precision of this measure- 
ment, therefore, will continue to improve over the course 
of Run II of the Tevatron. 
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FIG. 3: The confidence band in if* for a 95% C.L. is shown. 
Results from simulated experiments assuming a 955 pb~ 
dataset at different values of Tt are convoluted with a smear- 
ing function to account for systematic uncertainties. The fit- 
ted value from the data is indicated by an arrow. 
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